Magnetic Field scaling of Relaxation curves in Small Particle Systems 
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We study the effects of the magnetic field on the relaxation of the magnetization of small mon- 
odomain non-interacting particles with random orientations and distribution of anisotropy constants. 
Starting from a master equation, we build up an expression for the time dependence of the magne- 
tization which takes into account thermal activation only over barriers separating energy minima, 
which, in our model, can be computed exactly from analytical expressions. Numerical calculations 
of the relaxation curves for different distribution widths, and under different magnetic fields H and 
temperatures T, have been performed. We show how a rin(t/ro) scaling of the curves, at different 
T and for a given H, can be carried out after proper normalization of the data to the equilibrium 
magnetization. The resulting master curves are shown to be closely related to what we call effective 
energy barrier distributions, which, in our model, can be computed exactly from analytical expres- 
sions. The concept of effective distribution serves us as a basis for finding a scaling variable to scale 
relaxation curves at different H and a given T, thus showing that the field dependence of energy 
barriers can be also extracted from relaxation measurements. 



PACS Numbers: 75.1Q.Hk,75.4Q.Mg,75.50.Tt,75.60.Lr. 
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I. INTRODUCTION 



Time dependent phenomena in small-particle systems 
have been the subject of an increasing number of exper- 
iments because of their iiiterest as non-equilibrium phe- 
nomena in spipi systems jj for magnetic recording materi- 
als technologya and even as a possible way to prove ex- 
perimentally the existence of macroscopije |auantum tun- 
neling phenomena in magnetic materialsU^Q Whereas the 
basis of a theory of the magnetic aftM--|effect dates back 
from old studies on rock magnetismJ^iZI the interpreta- 
tion of several experimental results is still waiting for 
suitable theoretical models that capture the relevant fac- 
tors and parameters that can play a role in the explana- 
tion of these phenomena. One of the points that has not 
been completely clarified is the influence of a magnetic 
field in the relaxation of small-particle systems. 

Relaxation in zero field is usually analyzed in terms 
of parameters such aSr-jtiie so-called magnetic viscnfe 
ity S'js fluctuation fieldETlB and activation volumeJl3't^ 
which are susceptible|-tp_cnisinterpretations. In the last 
years, several authorsliTO have proposed an alternative 
method to analyze relaxation curves based on a T ln(t/ro) 
scaling of the relaxation data at different temperatures 
that avoids the above mentioned problems and gains in- 
sight on the microscopic details of the energjf barrier dis- 
tribution f{E) producing the relaxation.ll3ll2l In this con- 
text, the purpose of this article is to extend this kind of 
analysis to the case of relaxation in the presence of a 
magnetic field. We want to account for the experimen- 
tal studies on the relaxation of small-particle systems, 
which essentially measure the acquisition of magnetiza- 
tion of an initially demagnetised, punple under the ap- 
plication of a magnetic field Jl^'EirEd In this kind of ex- 



periments, the field modifies the energy barriers of the 
system that are responsible for the time variation of the 
magnetization, as well as the final state of equilibrium 
towards which the system relaxes. The fact that usu- 
ally the magnetic properties of the particles (anisotropy 
constants, easy-axis directions and volumes) are not uni- 
form in real samples, adds some difficulties to this anal- 
ysis because the effect of the magnetic field depencls-|(Mj. 
them in a complicated fashion. In a previous study jij'EZI 
we started to address some of these peculiarities, show- 
ing how experimental relaxation data must be treated 
in order to compare relaxation curves at different tem- 
peratures and fields making simple assumptions about 
the sample composition. Here, we will present the theo- 
retical background that supports this phenomenological 
approach, as well as detailed numerical calculations of 
the time dependence of the magnetization of a system 
of non-interacting randomly oriented small monodomain 
particles with uniaxial anisotropy and with a distribution 
of anisotropy constants. In a first approximation, we will 
neglect inter-particle interactions leaving for a future in- 
vestigation the effects of long-ranged dipolar interactions 
between the particles. 

The paper is organized as follows. In Sec. ||, we 
present the basic features of model show how the dis- 
tribution of energy barriers of the system is influenced b 
the application of a magnetic field with the help of the 
concept of effective energy barrier distribution. In Sec. 
[II we introduce the Two-State Approximation (TSA) for 
the calculation of the thermal dependence of the equilib- 
rium magnetization In Sec. IV, we derive the equation 



governing the time dependence of the magnetization from 
a master rate equation in the TSA. The results of numer- 
ical calculations based on the above mentioned equation 



1 



are presented in Sec.^. There, we present the Tln(t/To) 
scaUng of relaxation curves at a given magnetic field, dis- 
cussing its range of validity. We also study the possibility 
of a scaling at different fields and fixed temperature, and 
its applications. Finally in Sec. ^ we resume the main 
conclusions of the article. 



II. MODEL 

We consider an ensemble of randomly oriented non- 
interacting single-domain ferromagnetic particles of vol- 
ume V and magnetic moment M — MsVm with uniaxial 
anisotropy. To take into account the spread of particle 
volumes in real samples, we will assume that the parti- 
cles anisotropy constants K are distributed according to 
some function f{K). 

The energy of a particle is determined by the orien- 
tation of M with respect to the external magnetic field 
H and to the easy-axis direction n. Using the angular 
coordinates defined in Fig. it can be written as 

E ^ -— ^ -cos^{e) ~2hco8{e -iP) . (1) 
VK 

where we have defined the reduced field h = H/Hc and 
He = 2K/Mg as the critical field for an aligned particle. 
We have concentrated in the two dimensional case (M 
lying in the plane formed by H and n; ip — 0, since the 
energy maxima and minima can be calculated analyti- 
cally only in this case. In Fig. |l|, we show the variation 
of the energy with 9 for a typical case, defining in the 
same figure the notation for the energy barriers and ex- 
trema. 



A. Effective energy barrier distribution 

The magnetic field modifies the energy barriers of 
the system depending on the particle orientation and 
anisotropy value, and, consequently, changes the origi- 
nal energy barrier distribution.E3 Let be the energy 
barrier in zero field. Then, for a particle oriented at an 
angle "0, h mo difie s the barrier by a factor g{h, ip) in the 
following forn£a'E2l 

E,, = El g{h,i,) . (2) 

If f{ED is the energy barrier distribution in zero field, 
which has in fact the same functional dependence than 
the distribution of anisotropy constants f{K), then the 
distribution in the presence of a field is simply modified 
to 

fMh,E^,^) = f{El) (^^^) = f{El)/g{K^) , 

(3) 

which we will call effective energy barrier distribution. 



In order to understand the qualitative change of /off 
with h, we have numerically calculated fes{E\y) for a sys- 
tem of oriented particles with logarithmic-normal distri- 
bution of anisotropics 

f{K) = ^-ln-iK/K„)/2.^- 

for different widths a and Kq = 1, and several values 
of the magnetic field h. The calculation have been per- 
formed by making energy barriers histograms for a col- 
lection of 10 000 particles. The results are given in Fig. 
^ (upper panels). In all the cases, we observe the pro- 
gressive splitting of the original distribution in two 
subdistributions of high and low barriers as h increases 
from zero. The field tends to make deeper one of the 
minima, therefore increasing the two energy barriers for 
rotation of M out of the field direction, while the other 
two are reduced. In this way, the global effect of /i is a 
splitting of f{Eh) towards lower and higher values of Eh- 
As h attains the critical value for the particles with 
smaller a peak of zero or almost zero energy bar- 
riers starts to appear (see for example the curves for 
H — 0.5, 1.0 in the case a = 0.5); while most of the non- 
zero barriers are distributed according to a distribution 
identical to /{E^), but centered at higher energies. The 
higher the width of the distribution cr, the lower the h at 
which the lowest energy barriers start to be destroyed by 
the field. 

Finally, the combined effect of random orientations and 
f{K) has been considered. The results are shown in Fig. 
H (lower panels), where we can see that the features of 
the preceding case are still observed. Now, at high h, 
the distributions are smeared out by the disorder, and 
the minima becomes less pronounced due to the spread 
in particle orientations. 

In Sec. 0, we will discuss how these results affect the 
time dependence of magnetization in relaxation experi- 
ments. 



III. TWO-STATE APPROXIMATION 

The calculation of the equilibrium magnetization at 
non-zero T and finite K proceeds along the standard 
techniques of statistical mechanics. For particles oriented 
at an angle ^p, m(H, T) is simply given by the average of 
the projection of the magnetic moment of the particles 
onto the field direction over , all their possible orientations 
9. In our model, this isEZlfia 

m{H,T,^/j)=^ - [ dncoa9e-^^^^'l'\ (5) 
■2 Jn 

where is the solid angle and Z is the partition func- 
tion of the system. Here, the energy U {9, ip) appearing 
in the Boltzmann probability, has to be calculated from 
Eq. d), then 
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Ui0) = -a sin^ + ^ cos{0 - ^P) 



where the two dimensionless parameters 



(6) 



(7) 



have been introduced. 

At T such that the thermal energy fc^T is smaller 
than the relevant energy barriers of the system, typi- 
cally of the order of the anisotropy energy KV (a ^ 1), 
the main contribution to thermodynamic averages comes 
from states around the energy minima, since thermally 
activated jumps out of the stable directions of the mag- 
netization have extremely low probability to succeed. 
Therefore, as it will be useful for the numerical calcu- 
lations of the relaxation curves in Sec. 0, we will £Qftn 
sider the so-called Two-State Approximation (TSA).c3'E£l 
In this approximation, the continuum of states corre- 
sponding to all the possible orientations of m is truncated 
to the two local energy minima states. 

This will allow us to replace the integrations over mag- 
netization directions by sums over the two energy min- 
ima. If the particle has only one minimum, the two states 
considered in the calculation will be the minimum and 
the maximum of the energy function. For a system of 
randomly oriented particles and with a distribution of 
anisotropy constants f{K), Eq. (||) becomes in the TSA 

mTs{H,T)= dK dijfiK)mTs{K,ij), (8) 



where 



rriTS 



(9) 



i=l,2 



stands for the magnetization of an individual particle in 
the TSA, and /3= 1/A:bT. 

Eq. (^) has been numerically evaluated for a system 
of randomly oriented particles and several values of Kq 
and the results are displayed in Fig. ||. For the smallest 
Kq values, the curves present a small jump at a certain 
value of ^. This may seem unphysical but, in fact, this 
jump appears at an h equal to the critical field for the 
disappearance of one of the energy minima. In fact, when 
averaging over a distribution of anisotropics f{K) with 
Kq — 1 and tr — 0.5, this jump disappears. 

As expected, the TSA curves coincide with the re- 
sults obtained from the exact expression Eq. (^) for 
high enough Kq (compare the Kq = 10 case with the 
dashed-dotted line in Fig. On the other hand, at 

low enough Kq, the TSA reproduces the exact result for 
aligned particles, for which the magnetization curve re- 
duces to mxs = tanh(^), since the magnetization does 
not depend on a in this case, (compare the continuous 
line with the case Kq = 0.5) HI 




IV. RELAXATION CURVES IN THE PRESENCE 
OF A MAGNETIC FIELD 

Within the context of the Fokker-Planck equatiap^^ 
for M in the discrete orientation approximation,E3Ej we 
will assume that the relaxation of the magnetization due 
to thermal fluctuations can be modeled by a markovian 
stochastic process. Its dynamics can then be described by 
a master equation for Pi, the probability to find the mag- 
netization vector at time t in the equilibrium state i. Fur- 
thermore, we will assume that we are in the regime where 
the TSA is valid and, consequently, only transitions be- 
tween the two equilibrium directions of the magnetization 
given by the minima of the energy (^ will be considered. 
Moreover, in models considering continuous vaaables for 
the numerical evaluation of relaxation dynamicsEflLS, the 
elementary time step depends on T and H , giving rise to 
relaxation curves which are not directly comparable. In a 
recent work, Novak and ChantrelEll have faced the prob- 
lem of the quantification of the time step used in Monte 
Carlo simulation, giving a method to quantify the time 
step in real units. As an alternative, we propose a simple 
dynamical model that avoids this problem since, in the 
TSA, it can be solved analytically in terms of intrinsic 
parameters. 

Taking into account that the transitions between the 
two minima can take place cither by jumping over the 
barrier placed to the right or to the left of the initial state, 
the master equation governing the itipie dependence of 
the magnetization can be written asE3 



dt 

(k) 



EE{ 

k=l,2 j^i 



(10) 



where w'ly designates the transition rate for a jump from 
the state i to the state j separated by the maximum k 
(see Fig. |l|). The transition rates can be freely assigned 
as long as to fulfill the detailed balance condition.l23 It is 
a common choice to consider the Boltzmann probability 
with the energy difference between the two minima in 
the exponent. This choice, in spite of giving the correct 
thermodynamic averages in a Monte Carlo simulation, 
may not be appropriate to describe the dynamics of the 
system, since the energy barriers between the minima are 
not taken into account. 

For this reason, in the exponential of the Boltzmann 
probability, we have considered the energy difference be- 
tween the initial minimum i and the maximum k that 
separates it from the final state j 
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where Tq^ is the attempt frequency. It is a trivial mat- 
ter to prove that the following detailed balance equation 
holds 
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guaranteeing that thermal equilibrium is reached in the 
long time limit .l3 e = E^^^ — E^^^ is a measure of the 
asymmetry of the energy function. 

Taking into account the normalization condition Pi + 
P2 = 1, one can easily solve Eq. (jl^) for Pi and P2 as a 
function of time 



Pi{t) 



1 



o/3e 



-t/T\ 



1 



(13) 



The time-dependence of the system is thus characterized 
by an exponential function with a single relaxation time 
T that takes into account all possible probability fluxes 



= W = 



1 



(1 + e^') 



(14) 



As we see, r is dominated by the lowest energy barrier 
E^ , but with non- negligible pre- factors that take into ac- 
count the possibility of recrossing from the equilibrium 
to the metastable state and the two different possibili- 
ties of jumping. Notice that this two prefactors are often 
neglected in theoretical studies of the. dependence of the 
blocking temperature with the fielctl and Monte Carlo 
simulationsEJO. This is due to the fact that, usually, 
the possibility of jumping between minima by any of the 
two channels is not considered. However, at small non- 
zero fields (e > 0), and for particles oriented at -0 ^ 0, 
they can be equally relevant. This expression reduces to 
the usual one 



(15) 



when the energy function is symmetric (e = 0) and there 
is only one energy barrier, except for a factor 4 that can 
be absorbed in the definition of the prefactor tq. 

The time dependence of the magnetization of the par- 
ticle is then finally given by: 

m(t;X,V) = cos[0^„(7/.)]Pi(i) + cos[Ci„(V^)]P2(t) 

= niTsiK, tfj) + [mo - rhrsiK, V')]e-*/^(-^''^). (16) 

In this equation, mTs{K^^) is the equilibrium magne- 
tization in the TSA [Eq. (^], that has already been 
calculated in subsection III, and mo is the initial mag- 
netization. If we have an ensemble of randomly oriented 
particles and a distribution of anisotropy constants f{K), 
then the relaxation law of the magnetization is given by 



m{t) 



dKf{K) / dijj m{t;K,i>) 



(17) 



This will be the starting point for all the subsequent nu- 
merical calculations of the relaxation curves and mag- 
netic viscosity. 



V. NUMERICAL CALCULATIONS 

A. Relaxation curves: Tln(t/ro) scaling and 
normalization factors 

In this section, we present the results of numerical cal- 
culations of the magnetization decay based on Eq. ( [T^ ) 
for a system of particles with logarithmic-linear distribu- 
tion of anisotropy constants and random orientation of 
the easy-axis. For the sake of simplicity, we have assumed 
zero initial magnetization mo = 0, so particles have ini- 
tially their magnetic moments at random and evolve to- 
wards the equilibrium state mcq. In the following, we will 
use dimensionless reduced variables for temperature and 
time, defined ssT/Tq and t/ro, with Tq = Eo/kB and Eq 
the value of the energy at which f{K) is centered. 

We have assumed that the magnetic moment of each 
particle is independent of the volume, although, in fact, 
it can be proportional to it, but this effect can be eas- 
ily accounted by our model by simply changing f{K) 
to Kf{K) in all the expressions. For the case of 
a logarithmic-linear distribution, this change does not 
qualitjatiiiftljt modify the shape of the distribution. Other 
workgHjCj'Cil consider also a distribution of anisotropy 
fields He due to the spread of coercive fields in some real 
samples, but they study only relaxation rates at a fixed 
time. Here we have preferred to distribute K and the 
easy-axes directions, which has a similar effect, in order 
to separate as much as possible the effects of an applied 
magnetic field from other effects that may possibly lead 
to non-conclusive interpretation of the results. 

In Fig. |, we show the results of the numerical calcu- 
lations for a system with a = 0.5 for three different fields 
H — 0.1, 0.5, 1.0 and temperatures ranging from 0.02 to 
0.2. In the upper panels, we present the original relax- 
ations normalized to the equilibrium magnetization value 
as given by Eq. (|^). Normalization is essential ia order 
to compare relaxations at different temperaturealZI, espe- 
cially at low fields where the temperature dependence of 
the equilibrium magnetization is more pronounced. 

Our next goal is to investigate the possibility of scaling 
relaxation curves at different T in a given magnetic field 
with the scaling. variable Tlii{t/To), in the spirit of our 
previous workscjU. For this purpose, in the lower pan- 
els of Fig. ^, we show the relaxation curves of the upper 
panels as a function of the scaling variable Tln(i/To). Ac- 
cording to Ref. |l^, in absence of a magnetic field, scaling 
should be valid up to temperatures such that Te is of the 
order of a. Instead, we observe in Fig. ^ that, the higher 
the field, the better the scaling of the curves is in the long 
time region and the worse at short times. This observa- 
tion holds independently of the value of a, indicating that 
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it is a consequence of the application of a magnetic field. 
This can be understood with the help of the effe ctive en- 
ergy barrier distribution introduced in Sec. If A. As was 
shown in Fig. ||, h widens fcs{E) and shifts the lowest 
energy barriers towards the origin, giving rise to a sub- 
distribution of almost zero energy barriers that narrows 
with h, and, consequently, the requirements for Tln{t/To) 
scaling are worse fulfilled at small Tln(i/ro) values. On 
the contrary, as we will show in the next subsection, h 
broadens the high energy tail of energy barriers that con- 
tribute to the relaxation, f{E^^), improving the scaling 
requirements at large Thi{t/To) values. 



B. Scaling of relaxation curves at different magnetic 

fields 

Another interesting point is the possibility of finding 
an appropriate scaling variable to scale relaxation curves 
at different fields for a given T, in a way similar to the 
case of a fixed field and different temperatures, in which 
Tln(i/ro) is the appropriate scaling variable. In a first 
attempt, we will study the effect of ft- on a system with 
random anisotropy axes and the same K — 1. 



1. Randomly oriented particles, K = 1 



the barrier of a particle aligned with the field. In Fig. 
Pd, we have plotted the field dependence of all this quan- 
tities, together with the position of the maximum of the 
viscosity in energy units -Emax, and in Fig. |^ the value 
of the corresponding viscosity S'max- 

The reduction of tmax with h can be understood in 
terms of the progressive reduction of the energy barri- 
ers hy h. At /i = 0, the barriers are independent of the 
orientation of the particle and equal to 1, so that the 
maximum is placed at i?max = 1 and S'max = 1/e accord- 
ing to Eq. @. 

For h > 0.5 (the critical field for particles oriented 
at i/j = 45°), the lowest energy barriers start to be de- 
stroyed by h and consequently the relaxation rates peak 
at -Emax — with an increasing Smax value that increases 
as more particles loose their barriers. For /i > 1 all bar- 
riers have been destroyed and relaxations become field 
independent, with i?max = and S'max = 1/e. For fields 
up to h = 0.5, the variation of i?max and Smaxwith h 
can be used to scale the magnetic relaxation curves at 
constant T and different h. Therefore, although in this 
case the inflection points of the relaxation curves could be 
brought together by shifting them in the Tln(</To) axis 
in accordance with the variation of E'^'^, the full scal- 
ing cannot be accomplished because of the complicated 
variation of Smax (see Fig. ^). 



We have calculated the relaxation curves for this sys- 
tem at r = 0.05 and several values of the field. The 
obtained curves have been normalized to the equilibrium 
magnetization as given by Eq. (^). 

The effect of h on M (t) is better understood in terms 
of the logarithmic time derivative of M(t) 

which is the so-called magnetic viscosity S{t). As can be 
clearly seen from Fig. the viscosity curves at differ- 
ent h cannot be scaled neither by shifting them in the 
horizontal axis, nor by multiplicative factors, since the 
high and low field curves have different shapes. As soon 
as the field starts to destroy some of the energy barriers 
{h > 0.5), the qualitative form of the relaxation changes. 
This fact hinders, in principle, finding a field dependent 
scaling variable, valid in all the range of fields, in systems 
of non-aligned particles. 

Nevertheless, even though viscosity curves are qualita- 
tively different at different h, all of them present a well- 
defined maximum corresponding to the inflection point of 
the relaxation curves. This maximum appears at a time 
^max associated to an E'max = rin(tmax/''o), that de- 
creases with increasing h for a given temperature. This 
energy is approximately equal to the averaged lowest en- 
ergy barrier of the particles {E^^ in the notation of Sec. 
Q) and this value is closer to the lowest possible barrier 
(corresponding to particles oriented at i/" = 45°) than to 



2. Randomly oriented particles with f{K) 

In spite of the lack of scaling of the preceding case, 
in what follows, we will demonstrate that the inclusion 
of a distribution of K, always present in experimental 
systems, allows to scale the relaxation curves for a wide 
range of h. 

Let us consider a logarithmic-linear distribution of 
anisotropy constants of width a, Eq. (^. Low temper- 
atures relaxation rates corresponding to cr = 0.2, 0.5 are 
presented in Fig. ^. In this case, the qualitative shape 
of the viscosity curves is not distorted by h. It simply 
shifts the position of the maxima towards lower values 
of T\n{t/To) and narrows the width of the peaks, being 
these effects similar for both studied a. 

The position of the maximum relaxation rate still de- 
creases with increasing ft, following the decrease of the 
smallest energy barriers (see Fig. |^), which now have 
an almost linear dependency on ft. As in the preced- 
ing case, ii^max goes to zero when ft starts to destroy the 
lowest energy barriers. The difference is now that, due 
to the spread of the anisotropy constants, lower fields 
fto are needed to start extinguishing the lowest energy 
barriers (see Fig. |^), being this reduction greater, the 
greater a is, since the most probable anisotropy con- 
stant [_fi'max=^o exp(— cr^/2)] becomes smaller and, con- 
sequently, E'max drops to zero at smaller fto. This field 
corresponds to the one for which /eff(iJ^^) starts to de- 
velop a peak corresponding to zero energy barriers. As in 
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the preceding case, we have also tried to identify the vari- 
ation of E-aiax with the microscopic energy barriers of the 
system. As can be clearly seen in the dashed lines of the 
upper panels of Fig. the h dependence of £^max follows 
that of the lowest energy barriers for particles oriented 
ai ij: — 11 / 4, and with K = /^max- 

By looking in detail the low T relaxation curves {T — 
0.01 curves, analogue to the ones shown in Fig. ^ for 
T = 0.05), we have observed that two relaxation regimes 
can be distinguished. One presents a broad peak in S 
at relatively high energies (long times) with a maximum 
at an i?max which varies as the lowest energy barrier at 
^ = 45, this is clearly visible at the lowest h values even 
for the T — 0.05 case of Fig. The other regime presents 
a peak around E — that starts to develop as soon as h 
breaks the lowest energy barriers. What happens is that 
the first peak shifts towards lower energies with h, at the 
same time that the relative contribution of the second 
peak increases, the global effect being that, at a certain 
h, the contribution of the first peak has been swallowed 
by the second because at high h and low T, relaxation is 
driven by almost zero energy barriers. 

To clarify this point, we show in Fig. || S{t)/T for 
three different temperatures and magnetic fields H = 
0.1,0.5, 1.0 for a narrow {a ~ 0.2) and a wide {a — 0.5) 
distribution. The effective distribution of lowest energy 
barriers /cff(i5^^), already calculated in Fig. ^, is also 
plotted as a continuous line. We observe that for a nar- 
row distribution, at low enough T, S{t)/T coincides with 
fcs{E'^^) independently of h, demonstrating that only 
the lowest energy barriers of the system contribute to 
the relaxation. 

Finally, let us also notice that, at difference with the 
preceding case, S'max becomes almost constant below ho 
(lower panels in Fig. ^ and low enough T, so that now 
the relaxation curves at different h and fixed T may be 
brought to a single curve by shifting them along the 
Tln(i/ro) axis in accordance to the i?max variation. The 
resulting curves are displayed in Fig. || for a — 0.20.5. 
They are the equivalent of the master curves of Fig. Wfor 
a fixed h and different T. Now the appropriate scaling 
variable is 

£;sca = Tln[t/W(M] , (19) 

which generalizes the scaling at fixed T. This new scal- 
ing is valid for fields lower than ho, the field at which the 
lowest barriers start to be destroyed and above which 
the relaxation becomes dominated by almost zero energy 
barriers. Thus, as already discussed in the previous para- 
graphs, the wider a, the smaller the h range for the va- 
lidity of field scaling. 

VI. CONCLUSIONS 

We have proposed a model for the relaxation of small 
particles systems under a magnetic field which can be 



solved analytically and which allows to study the effect of 
the magnetic field on the energy barrier distribution. In 
particular, we have shown that the original /(-Bb) is split- 
ted into two subdistributions which evolve towards higher 
and lower energy values, respectively, as h increases. It 
is precisely the subdistribution of lowest energy barriers, 
the one that completely dominates the relaxation as it is 
evidenced by its coincidence with the relaxation rate at 
low T. 

For fields smaller than the critical values for the small- 
est barriers, the relaxation curves at different h and fixed 
T can be collapsed into a single curve, in a similar way 
than rin(i/To) scaling for curves at fixed h. Whereas 
the latter allows to extract the baccier distribution by 
differentiation of the master curvdiZi, the shifts in the 
Tln(i/ro) axis necessary to produce field scaling, give the 
field dependence of the mean relaxing barriers, a micro- 
scopic information which cannot easily be inferred from 
other methodsLJ. 
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FIG. 1. Energy function E{8,tp) as a function of the angle 
between the magnetization vector m and the magnetic field h, 
for m in the plane of the easy-axis <^ = 0, as given by Eq. (|l[). 
The plot is for a particle whose easy-axis n forms an angle 
4> = 30° with h, and H — 0.3. We have used the following 
notation to designate the extrema of the energy: B^in ^-nd 
^max refer to the extrema closer to the field direction while 
^min Bind ^max refer to those further from the direction of the 
field. The four possible energy barriers between them are 
El'' = E{9lnax) — ^(^min)- Inset: Schematic representation of 
the quantities involved in the definition of the system. The 
easy-axis of the particles n are in the x-z plane forming an 
angle ip with the magnetic field H, which points along the z 
axis. and ip are the spherical angles of the magnetization 
vector M. 
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FIG. 2. Upper panels: effective energy barrier distribu- 
tions for aligned particles with a lognormal distribution of 
anisotropy constants of a = 0.2 (continuous lines), 0.5 
(dashed lines), 0.8 (dot-dashed lines) for values of H as indi- 
cated in the figures. Lower panels: same as upper panels but 
for particles with random orientations of anisotropy axes. 




FIG. 3. Magnetization curves as a function of the dime- 
sionless Zeeman energy ^ = ^HV/knT in the TS approxi- 
mation. Symbols stand for randomly oriented particles with 
Ko = 0.5, 1.0, 10.0 (from the uppermost curve). The case 
A'o = 10 is compared to the exact result given by Eq. 
(^) (dash-dotted line). The case Ko = 1 is compared with 
a system of randomly oriented particles with f{K), Ko = 1 
and a = 0.5 (long-dashed line) The result for aligned particles 
is displayed as a continuous line, for which rriTS ~ tanh(^) 
(rriTS is independent of a in this case). 
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FIG. 4. Relaxation curves for an ensemble of particles with randomly oriented anisotropy axes and a logarithmic-normal 
distribution of anisotropies f{K) of width a — 0.5 and Ko — 1 calculated by numerical integration of Eq. (p^). The initial 
magnetization has been set to Mo — 0. Reduced temperatures T/To, starting from the lowermost curve, range from 0.01 to 0.1 
with 0.01 increments, and from 0.1 to 0.2 with 0.02 increments. The applied fields are H — 0.1, 0.5, 1.0 as indicated. The upper 
panels show the original curves normalized to the equilibrium magnetization mTs(T) given by Eq. (^). In the lower panels, 
the same curves have been plotted as a function of the scaling variable rin(t/ro). 
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FIG. 5. (a) Low temperature (T — 0.05) viscosity curves 
for a system of randomly oriented particles with the same 
anisotropy constant Ko = 1- The curves have been normal- 
ized to the equilibrium magnetization Meq (T) and correspond 
to magnetic fields H = 0.1, 0.2, 1.0, 1.2, 1.4, 1.6, 1.8, 2.0 
increasing from right to left. (b) Field dependence of 
the time corresponding to the maximum relaxation rate, 
T \n{tmax / To) , Bs derived from the viscosity curves in panel 
(a) (triangles). The field dependence of the mean lowest en- 
orgy barrier e!^,j^^ (diamonds), lowest energy barrier for par- 
ticles oriented at a = 45°, i5'-^'^-'(7r/4) (squares) and a = 
(circles), E^^''^\Q°) are also shown for comparison, (c) Field 
dependence of the maximum relaxation rate Smax- 




0.0 0.5 
Tln(t/T:„) 
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FIG. 6. (a) Low temperature (T = 0.05) viscosity curves 
for a system of particles with random orientations and loga- 
rithmic-normal distribution of anisotropies with (a) a = 0.2 
and (b) a = 0.5. The curves have been normalized to the 
equilibrium magnetization and correspond to magnetic fields 
H = 0.1 to 1.0 in 0.1 steps and H = 1.2 to 2.0 in 0.2 steps 
starting from the right. 
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FIG. 7. Upper panels: Field dependence of the energy cor- 
responding to the maximum relaxation rate, T\n{tmax/To), 
as derived from the viscosity curves in Fig. ^ for temper- 
atures T — 0.01 (circles), 0.05 (squares) , 0.1 (diamonds). 
Lower panels: Field dependence of the maximum relaxation 
rate Smax for the same curves and temperatures. Left column 
is for a — 0.2 and the right one for a = 0.5. 
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FIG. 8. Relaxation rates as a function of the scaling vari- 
able Tln{t/To) for different temperatures [T = 0.01 (dotted 
line), 0.05 (dashed line) , 0.1 (long-dashed line)], a — 0.2 
and three magnetic fields (a) H = 0.1, (b) H — 0.5, and 
(c) H = 0.8. The curves tend to the effective distribution of 
lower energy barriers fcs{E'^^), shown as a continuous line, 
as T decreases. 



11 



0.0 
0.: 



.50 0.75 1.00 1.25 1., 



.50 



1.0 




0.00 0.50 1.00 1.50 2.00 

Tln[t/t„,,,/h)] 

FIG. 9. Normalized relaxation curves as a function of the 
scaling variable Tln[{t/tmax{h)] for T = 0.05) obtained from 
Fig. y by shifting the curves in the horizontal axis with the 
position of the maximum relaxation rate (upper panels in Fig. 
0). (a) a = 0.2 and H = 0.1, 0.2, 0.3, 0.4, 0.5; (b) a = 0.5 and 
H = 0.1, 0.2, 0.3, 0.4 (starting from the uppermost curve) 
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